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ABSTRACT 

We investigate the scaling relations between the X-ray and the thermal Sunyaev-Zel'dovich Effect 
(SZE) properties of clusters of galaxies, using data taken during 2007 by the Y.T. Lee Array for 
Microwave Background Anisotropy (AMiBA) at 94 GHz for the six clusters A1689, A1995, A2142, 
A2163, A2261, and A2390. The scaling relations relate the integrated Compton-j/ parameter I2500 to 
the X-ray derived gas temperature T c , total mass M2500, and bolometric luminosity Lx within r 2 5oo- 
Our results for the power-law index and normalization are both consistent with the self-similar model 
and other studies in the literature except for the Y 2 5oo _ Ljf relation, for which a physical explanation 
is given though further investigation may be still needed. Our results not only provide confidence for 
the AMiBA project but also support our understanding of galaxy clusters. 



Subject headings: cosmic microwave background — cosmology: observations 
general — X-rays: galaxies: clusters 
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1. INTRODUCTION 

The Sunyaev-Zel'dovich Effect (SZE) is a powerful 
tool that can potentially answer long-standing questions 
about the large-scale distribution of matter. The SZE is a 
spectral distortion of the cosmic microwave background 
(CMB), induced when a fraction of CMB photons are 
scattered by hot electrons in the core s of massive galaxy 
clusters (jSunvaev fc Zel'dovichl 11970ft . The redshift in- 
dependence of the SZE enables the direct detection of 
distant clusters without the (1 + z) 4 surface brightness 
dimming that limits other techniques, including X-ray 
observations. Clusters studied via the SZE are therefore 
effective cosmological probes. Studying their properties 
in detail will lead to heightened understanding of the 
mass power spectrum, and should provide improved con- 
straints on cosmological parameters. 

In the simplest scenario, where gravity is assumed to 
be the only influence on the formation of galaxy clus- 
ters, a simple 'self-similar' model can b e used to re - 
late the physical properties of clusters (jKaiserl fl986L 
Assuming spherical collapse of the dark matter (DM) 
halo, and hydrostatic equilibrium of gas in the DM 
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gravitational potential, one can derive power-law scal- 
ing relations between various X-ray and SZE quantities, 
e.g. luminosity and temperature, gas mass and temper- 
ature, total mass and luminosity, entropy and tempera- 
ture, and Compton-y parameter and temperature. Ex- 
istence of th ese relations in observ a tions c an be seen in, 
for ex a mple. iMushotzkv fc Scharil (I1997T): iBenson et al.l 
(|200l : iBonamente et all ([20081) : iMorandi et all (|2007l) . 
These relations are also found in numerical simula - 
tions, e.g., between X-ray quantities dNagai et al.ll2007ft . 
SZE flux an d total mass or gas mass (jMotl et al.l 120051 : 
INaeail 120061 ) , and in tegrated Compton-ji/ and temper- 
ature or luminosity (|da Silva et al.l |2004[ ) . Deviations 
from the scaling relations should reveal the importance 
of non-gra vitational processes for the formation of clus- 
ters (e . g. lAllen fc Fabianl [19981 : [McCarthy et al.l I2002L 
I2003al). or constrain the ma ss distributions of clusters 
(Rciprich & Bohringcr 2002). Furthermore, the scaling 
relations can be used as an utility to extract important 
quantities and evolution behaviors for remote clusters us- 
ing SZE observables alone. 

The Y.T. Lee Array for Mic rowave Backgr ound 
Anisotropy (AMiBA) experiment (|Ho et alj 120091 ) ob- 
served and detected the SZEs of six massive Abell clus- 
ters in the rang e 0.09 < z < 0.32 during the year 2007 
(|Wu e t al. 2009). AMiBA is a coplanar interferometer 
that during 2007 operated at 94 GHz with seven 0.6-m 
antennas in a hexagonal close-packed configuration, giv- 
ing a synthesized resolution of about 6'. The array has 
a sensitivity of 63 mJy/hr for on-source integr ation, and 
an overall efficiency of 0.36 (jLin et al.l [2009) . Details 
for the transformation of th e raw dat a into calibrated 
visibilities are presented in IWu et all (12009ft. and the 
checks on data integrity are described in iNishioka et al.l 
(2009). At our observing frequency the SZE signal is 
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an intensity decrement in the CMB. We fit the central 
(peak) decrem ent in the AMiBA visibilities using isother- 
mal /3- models (|Cavaliere fc Fusco-Femianolll976l ). taking 
account of contaminat ion from the p rimary CMB and 
foreground emissions (iLiu et al ir2009D Other compan- 
ion pap ers include iChen et al.l (|2009h and iKoch et al.l 
(2009a]), wher e the technical aspect s of the instruments 
are described, lUmetsu et~aTl l|2009|) . where the AMiBA 
SZE data is combined with weak lensing data from 
Subaru t o analyze t h e dist ributions of mass and hot 
baryons, IKoch et al.l (|2009bl ). where the Hubble con- 
stan t is estimated from AMiBA SZE and X-ray data, 
and iMolnar et al.l (|2009l ), which discusses the feasibil- 
ity of further constraining the intra-clu ster gas model 
using AMiBA upgraded to 13 antennas (|Ho et al.l I2009L 
AMiBAl3). The consistency of our results with other ob- 
servations and theoretical expectations will validate not 
only the performance and capability of instruments, but 
also the analysis methodology. Since AMiBA is one of 
few leading SZE instruments operating at 3-mm wave- 
length, we anticipate that it will fill an important role by 
providing 3-mm SZE data for spectral studies. 

In this article we address the scaling relations between 
the integrated SZE Compton-y parameter obtained by 
AMiBA and X-ray gas temperature, X-ray luminosity, 
and total cluster mass derived from the literature. In 
Section [2] we discuss the cluster gas models and cluster 
parameters derived from the X-ray data. In Section [3] we 
calculate the integrated Compton-y parameter for each 
of the six clusters. In Section[4]we investigate the scaling 
relations including the consideration for errors. We fur- 
ther discuss our results in Section[5]and draw conclusions 
in Section [6l 

2. CLUSTER PROPERTIES FROM X-RAY DATA 

As the u-v coverage is incomplete for a single interfero- 
mctric SZE experiment, we can not measure the accurate 
profile of a cluster or its central intensity. Therefore we 
have chosen to assume a cluster model, and thus a flux- 
density profile, so that a corresponding template in u-v 
space can be fitted to the observed visibilities in order 
to estimate the underlying model parameters including 
the central SZE intensity, A Jo. We apply the spherical 
isothermal /3-model in our X-ray and SZE analysis. The 
cluster gas density distribution is of the form 

/ r 2 X -3/3/2 

n c (r) = n c0 f 1 + J , (1) 

where n e o is the central number density of electrons, r 
is the radius from the cluster center, r c the core radius 
and B is a structure index. Due to the limited resolution 
of AMiBA in its 7-element closed-packed configuration, 
we cannot obtain good estimates for some of the model 
parameters from our SZE data alone. Therefore we have 
taken the X-ray derived values for 8 and 9 C from the lit- 
erature, where 9 C — r c /Dx and Da is the angular diam- 
eter distance. Throughout this paper, we assume a flat 
ACDM universe with H = 70 km s _1 Mpc _1 , fl M = 0.3, 
and J7a = 0.7. 

To relate with the SZE Compton-y parameter, we also 
need to borrow the cluster gas temperature T e , the total 
mass M2500, and the bolometric luminosity Lx derived 
from the X-ray data. M2500 refers to the total mass in a 



cluster central region out to r = f25oo , defined as the ra- 
dius of mean overdensity 2500 x p c where p c is the critical 
density at redshift z. Given B, r c and T e for a cluster, 
we can compute first r25oo and the n Afeon through th e 
total mass equation of the /3-model (jGrego et al.ll200lh 

M 2500 = 25 00 Y r 250oPc = ^ ^ - ^- , (2) 

where /i = 0.6 is the mean molecular weight in units of 
ttih for an assumed near-solar metallicity of the intra- 
cluster medium. 

We considered two sets of X-ray derived parameters. 
The first set is mainly based on the Chandra data, 
and this leads to our main results. The second set is 
mainly derived from ROSAT images or a combination 
of ROSAT data and ASCA spectral measurements (the 
ASCA/ROSAT parameters in what follows). Because 
these data are generally of lower accuracy, we include 
them only for comparison. 

2.1. Chandra 

To deal with the complicated non-gravitational physics 
in cluster cores, including radiative cooling and feedback 
mechanisms, and the transient boosting of surface bright- 
ness and spectral temperature during merging events, the 
parameters of the Chandra set were derived by fitting an 
isothermal /3-model to the X-ray data with the central 
100 kpc excised. The most recent and curren tly most 
extensive studies of H n dBonamente et al.ll2006|) and the 
gas mass fraction / gas ()LaRoaue et al.ll2006D adopted this 
100-kpc cut model in their analysis, and claim that a cut 
at 100 kpc is large enough to exclude the cooling region 
in cool-core clusters while retaining sufficient photons for 
modeling. This model was also used in recent studies of 
scaling relations based on X-ray and SZE o bservations 
(|Morandi et alj|2007t iBonamente et al.ll2008f) . 

Table [T] summarizes the parameters from Chandra ob- 
servations, and the values of r25oo and M250Q derived 
from them. Note that C/iandra-based values of 8 and 
9 C for A2142 are unavailable in the literature, and so we 
adopted values taken from the ASCA/ROSAT set, which 
were not determined by fitting the 100-kpc cut model. 
The gas temperature for A2142 is Chandra-based, as 
given by iMarkevitch et al.l (|2000l ). The temperature fit 
allowed a cooling component to be present, but was based 
on the overall X-ray spectrum of A2142, rather than dis- 
carding photons extracted from the central 100 kpc re- 
gion. Nevertheless, if A2142 is excluded from the sample 
for this set of parameters, it has only a minor effect on 
the scaling relations (less than a 5% change for the power 
index or the normalization; see Sec.[4|). 

The valu e s of 8 and 6 C for A2390 are taken from 
lAllen et all (|2001[ ) who fit the X-ray surface brightness 
profile to an isothermal /3-model between radii 80 and 
900 kpc. As the authors remark, however, an isothermal 
/3-model ignoring the central region associated with the 
possible cooling flow cannot describe the mass distribu- 
tion well, since there is a 'break' in the surface brightness 
profile at r ~ 500 kpc. A better fit can be obtained using 
a simple broken pow er-law model, or assuming an NFW 
(|Navarro et al.l 119971 ) potential with the assumption of 
gas isothermality relaxed. 
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2.2. ASCA/ROSAT 

Parameters derived from ASCA and ROSAT are sum- 
marized in Table [2] The gas temperatures and the bolo- 
metric luminosities of our clusters, except A 199 5 and 
A2163, are compiled bv lAllen fc Fabian! (| 1998ft and|Allen| 
(2000), where the X-ray spectra were fitted by using a 
model with an isothermal plasma in collisional equilib- 
rium, including an additional component explicitly to ac- 
count for cooling flows (Model C). For A2163, which is 
not a cooling-core cluster, we take the values from the 
same papers, but without the additional cooling compo- 
nent (Model A). For A1995, which i s absent from thes e 
papers, we use the value of T e from iPatel et alJ ((2000), 
who detected no excess in the X-ray surface brightness bi- 
ased from a cooling flow in the cluster center. However, 
A1995 has recently been classified as a cooling cluster, 
according to the criterion that the cooling time in the 
central inner re gion is less than the Hubble time at the 
cluster redshift (jMorandi et al.ll2007ft . 

All values in Tables [Hand [2] are presented at the 68.3% 
confidence level. Errors are obtained by propagating 
the errors in the input parameters from the literature 
through a Monte-Carlo process. 

3. CLUSTER PROPERTIES FROM SZE 

In AMiBA targeted observations at 94 GHz, the sky 
signal is dominated by the thermal SZE. The amplitude 
of such signals is proportional to the Compton-y pa- 
rameter, y — <jt I (m c c 2 ^j J k^T e {Y]n e (l) dl where ot is 
the Thomson scattering cross section, k-B,T e (l)n e {l) is the 
electron pressure, and the integral is taken along the line 
of sight. The Compton-j/ parameter can be interpreted 
as a measure of Comptonization integrated through a 
cluster. In terms of a change in intensity, the thermal 
SZE observed at frequency v can be represented by a 
decrement 



ZE 



g (x, r ) • Icmb , 
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where x = hv j (fc B TcMB), T C mb 
(jMather et al.lfl999ft . and J C mb = 2hv 3 c~ 2 (e x - 1) 
is the CM B intensity. The factor q (x, T e ) can be ex- 
pressed as (IBonamente et all 12551 iMorandi et all 120071 ; 
lUdomprasert et al.ll2004l ) 
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where S rP ^(x,T e ) is a s mall relativistic correction 
(jChallinor fc LasenbvllilM ) 
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F = xcoth(a;/2), and G = xj sinh(x/2). The relativistic 
correction is about 6% for v = 94 GHz and T e — 10 keV, 
which is a typical temperature for our SZE clusters. 

Given a gas density profile n e (r) we can determine 
the distribution of A/sze on the plane of the sky. 
For an isothermal /3-model, the projected SZE decre- 
ment distribution has a simple analytical form (e.g. 



lUdomprasert et al.ll2551 



a2 \ (l~3/3)/2 

A/sze(0) = A/ o (1+02 



(6) 



where 9 and C are the angular equivalents of r and r c 
respectively, and AIq is the central SZE intensity decre- 
ment. Because the SZE clusters are not well resolved by 
AMiBA, we cannot get a good estimate of AIo, (3, and 9 C 
simultaneously from our data alone. Instead, we adopt 
the X-ray derived values for j3 and 6 C f rom Chandra or 
ASCA/ROSAT, and then estimate AI (|Liu et al.ll200ll 
by fitting the /3- model to the SZE visibilities obtained in 
IWu et all (|25oa) . 

For the two different sets of X-ray parameters we ac- 
cordingly obtain two sets of AIq values. For the the 
Chandra set with a 100 kpc-cut model we choose to 
fit the entire SZE data, whi le using the X- r ay pa ram- 
eters from the same mo del lLaRoque et al.l (|2006ft and 
iBonamente et all ([2006) already remarked that there is 
no simple way to mask the central 100 kpc from the in- 
terferometric SZE data because these data are in the 
u-v space. Nevertheless, our approach should be valid 
because the limited resolution of the current AMiBA 
is insensitive to the details of the cluster core. More- 
over, since the SZE probes the integrated gas pressure, 
which is linear in n c , the parameters derived from the 
SZE data should be less dependent on the core proper- 
ties than parameters derived from X-ray observations, 
where the X-ray surface brightness oc n 2 . Table [3] sum- 
marizes the resulting estimated values of AJ based on 
Chandra. We note that the effects of foregrounds such as 
radio source contamination, Galactic emission, and con- 
fusion from p rimary CMB fl uctuations have been taken 
into account (|Liu et al.H2009ft . 

In addition to the intensity decrement, the thermal 
SZE also can be expressed in terms of a change of 
the thermodynamic temperature of th e CMB, ATsze = 
y-q (x , T c ) (e x — 1) (xe x )~ 1 -Tomb (e.g. . IBonamente et all 
2008). Thus for a cluster observed at a given frequency, 
the A/sze (Eq. ©) and ATsze are equivalent measures 
of the Compton-y parameter. In Tabled we include the 
values of the central thermodynamic temperature decre- 
ment ATo that correspond to the AIq based on Chandra. 

To obtain an overall measure of the thermal en- 
ergy content in a cluster, we computed the integrated 
Compton-y parameter I2500 1 which is the Compton-?/ in- 
tegrated from its center out to the projected radius r25ocb 
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where is the solid angle of the integrated patch and 
fi 2 5oo is the total value covered within radius T25oo • The 
integrated Compton-y parameter has been shown to be a 
more robust quantity than the central value of Compton- 
y for observational tests, because it is less dependent 
on the model of gas d istribution used for the analy- 
sis (|Benson et al.l [2004). In addition, integrating the 
Compton- j/ out to a large projected radius diminishes 
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TABLE 1 

Cluster parameters of the Chandra set 
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A1689 

A1995 
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A2261 

A2390 


0.183 
0.322 
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0.224 
0.233 


0.686±°;°* 

003+0 021 
u.»zo_ 023 

74+ 01 
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n co+0.058 
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z - D '-0.32 
c 7 +0.61 
°'°'-0.58 


3.15 ±0.09 
1.51 ±0.05 

4.80 ±0.05 
2.02 ±0.07 
4.66 ± 0.05 


1, 4, 4 

1, 4, 4 

2, 5, - 
1, 4, 4 
1, 4, 4 

3, 4, 4 



Sanderson & Ponman (2003); Lancaster ct al. (2005). 
( 20071) . (51 IMarkeviteh et al.l l|2000f l . 



(3) lAllen et afl 



Ref erences. — (1) [Bonamcntc et al. (200 6j). (2) 
||2001|) with a 10% error assumed. (4) Morandi ct al 
a The rodshifts z are from IBonamente et all ( 1200611 except those for A2142 & A2390 which are given by I Alien! (120001) . 
b We take the values of and 6 C used in the ASCA/ROSAT set (Table [5J for A2142 since they are not available in the 
Chandra-based literature. 

c The emission-weighted temperatures and the bolomctric luminosities are extracted in a region of radius r between 100 
kpc and BOO- 



TABLE 2 

Cluster parameters of the ASCA/ROSAT set 
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14.7 
5.83 
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1, 4, 4 

1, 5, - 

2, 4, 4 
1, 4, 4 
1, 4, 4 

3, 4, 4 



Ref erences. — (1) IReese et al.l (f20o3) . ( 21 [Sanderson fc Ponmanl (1200311 : [Lancaster et al.l (1200511 . (31 IBohringer et al.l 
(1998) with a 10% error assumed. (4)lAllen (200011. (5) IPatel et alj ( 120001) 



(though does not completely remove) effects resulting 
from the presence of str ong entropy features in the cen- 
tral regions of clusters (jMcCarthv et al.ll2003af) . Table 
[3] summarizes our derived values of Y2500, adopting the 
parameters based on Chandra. In Sectional the Y2500 de- 
rived from both X-ray parameter sets will be considered 
for its scaling relationship with T e , M2500 and Lx- 

Although using X-ray data to determine the shapes of 
cluster SZE profiles is a common strategy in SZE analy- 
sis, it has been shown that this will bias the results of fit- 
ted parameters due to the assumption of i sothermality of 
a /3-m odel (e.g. IKomatsu &: Seliakll200U Hallma n et al.l 
120071) . In Section [5] we will further discuss this issue, 
and ap ply a simple correction to our results based on the 
work of lHallman et alj (|2007D . 

AMiBA is one of the first instruments to pro- 
vide 3-mm SZ detections of the cluster targets, ex- 
panding our knowledge of the SZE spectra for clus- 
ters. Table [3] compares our results for Y2500 at 
94 GHz with results at ot her frequencies: the 
BIMA/OVRO r esults at 30 GHz ([McCarthy et al.ll2003bl ; 
iMorandi et"aH 120071) . and th e SuZIE II results at 
145 GHz (jBenson et al.l I2004D . We have converted 
the BIMA/OVRO valu es of integrated Compton-y pa- 
rameter, 2/2500 (|Morandi et al.l 120071 ). to our Y2500 
using F 250 o = 2/2500 (F - 4) xe x / [(e x - 1) I g(x, T e )] 
where Iq = 2 (fcB?CMB) 3 {hc)~ 2 , and x, F, and 
g(x,T e ) are defined in Eq. ([4]) with frequency v = 



30 GHz. For SuZIE II, F 2 500 is obtained from y 2500 = 
S (r 25 oo) / [/cmb g(x, To) }, where S (r 2 Rno) i s the inte- 
grated SZ flux defined in lBenson et al.1 (|2004| ). and Jcmb 
and g(x, T e ) are calculated at v = 30 GHz. All three sets 
of results are based on reconstruction of the gas profile of 
the clusters using an isothermal /3-model, and all include 
the relativistic correction in their estimates of ^2500- Our 
results are consistent with those from BIMA/OVRO ex- 
cept for A1995, and can be seen to be generally lower 
than those from SuZIE II. 

Since Y2500 is relatively insensitive to cluster mor- 
phology and the thermal structure of the intraclustcr 
medium, the difference between the results for A1995 
could mostly come from the different SZ techniques. The 
AMiBA data led to a signal-to-noise ratio of about 6, 
based on an integration time of about 5.5 hours. A1995 
is a relatively low-mass, cool, cluster and hence produces 
a smaller SZE than the other clusters in our sample. This 
causes the statistical significance of the different values 
for y 2 5oo to be low, and that differing residual contam- 
ination by point sources and primary CMB fluctuations 
(|Liu et alj 120091 ) could be an important factor. Further 
multi-frequency studies of A1995 are needed to improve 
the result for the SZE of this cluster. 

4. SCALING RELATIONS 
4.1. Theoretical Expectations 
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TABLE 3 

Parameters of AMiBA clusters derived from SZE observations 



Cluster 


AJ 

(10 5 Jy/sr) a 


AT 
(mK) 


Y 2 

BIMA/OVRO 


500 (10- 10 sr) 
AMiBA 


SuZIE II 


A1689 
A1995 
A2142 
A2163 
A2261 
A2390 


-2.36 ± 0.71 
-3.19 ± 1.23 
-2.09 ±0.36 
-3.64 ±0.61 
-2.59 ± 0.90 
-2.85 ± 0.77 


-0.40 ±0.12 
-0.54 ±0.21 
-0.35 ± 0.06 
-0.62 ±0.10 
-0.44 ±0.15 
-0.48 ±0.13 


2.17 ±0.14 
0.71 ± 0.06 

5.53 ± 0.41 
1.51 ± 0.18 


2.82 ±0.86 
1.49 ±0.58 
13.44 ± 2.40 
6.61 ± 1.38 
1.72 ±0.64 
3.12 ±0.98 


465 +0.61 

5 50+°- 76 

4 - 46 -0.94 

3 6Q+ ' 56 
o.oy_ 57 



Note. — The integrated Compton parameters Y2500 measur ed by AMiBA (94 GHz) 
are compared with results from BIMA/OVRO (30 GHz: IMcCarthv et "aU I2003bt 
Morandi c t al.l2007f) and SuZIE II (145 GHz; deduced from lBenson et al.ll2004fl . The cen- 
tral SZE intensity AIq, its corresponding thermodynamic temperature decrement A To, 
and the AMiBA Y2500 were derived using the Chandra-based parameters. Isothermal 
/3-modcls arc used in all three sets of observation to reconstruct the gas profile of clus- 
ters and derive Ybsoo- The relativistic correction 5 re i (x,T c ) in Eq. |5j is also taken into 
account in all three cases. Errors are given at the 68.3% confidence level. 
a Central SZ intensities arc given by Liu ct al. (2009). 



In the context of the self-similar model, if assuming 
hydrostatic equilibrium and an isothermal distribution 
of baryons in the spherically-collapsed DM halo, it can 
be shown that there are simple power-law scaling rela- 
tions between the SZE and X-ray quantities. Specifically 
there are simple relations between the integrated Comp- 
tonization and the gas temperature T , the cluster total 
mass M to t and the bolomctric X-ray luminosity Lx, 

YD\ ocT? /2 E(z)-\ (8) 

YDl^MU?E{ Z yl\ (9) 

YDlotL^Eiz)- 9 / 4 , (10) 

where E 2 (z) = Q M (1 + zf + Q A + ^k(l + z) 2 
(jMorandi et alJ 120071 ) . We note that these scaling re- 
lations assume that the fraction of the cluster mass 
presen t as gas, / gas , is a constant. iBonamente et all 
(2008) found no significant scatter of f gas in their re- 
sults. Nevertheless, recent X-ray work in observations 
(e.g. IVikhUniril 120061 ) and simulations (jKravtsov et alJ 
2005) suggest that some variation m ay be expected. 

Following standard method (e.g. [Press etail f2002Tl . 
we perform a linear least-squares fitting in log 10 space, 
log 10 (y) = A + B\og 10 (x), taking account of errors in 
both x and y, to estimate the normalization A and power 
law index B of each scaling relation. The \ 2 statistic is 
defined as 

(a Vi log 10 (e)/yi) 2 + (Ba Xt log 10 (e)/xi) 2 

(11) 

and is minimized as in IBenson et ail (|2004 Eq. (13)). 
a Vi and a Xi for the sample points are obtained from the 
upper and lower uncertainties around the best-fit values 
as a = (cr + + er~) /2. The number of degrees of freedom 
(d.o.f.) is N — 2 with N equal to the total number of clus- 
ters in the sample, la errors in A and B are determined 
by projecting the A% 2 = 1 contour on each coordinate 
axis. 

4.2. Derived Observational Results 

The results of fitting log 10 (y) = A+B logjr^x) for each 
scaling relation are summarized in Table [4] Figures [TJ 



[21 and [3] show our sample of six clusters and the best- 
fitting scaling relations. Each figure shows the scaling 
results based on the Chandra and the ASCA/ROSAT 
parameters, for comparison. Five of the six clusters in 
our sample are cooling-core (CC) clusters; the exception 
is cluster A21 63, which is of non - cooling core (NCC) type 
(iMvers et alJ ll^jAllenl [2001 IMcCarthv etail [2003b; 
iMorandi et al.ll2007l ). " 

4.2.1. The Y2500 -T e relation 

Our results for the power law index B from both sets 
of X-ray parameters agree with the self-similar model 
B = 2.5 at the la level. They are also consistent 
with the values of B = 2.37 ± 0.23 from BIMA/OVRO 
( Bonamente "eTal1l2 008) . B = 2.21 ± 0.41 from SuZIE II 
( Benson et all 12004 . and B = 2.64 ± 0.28 (CC+NCC 
sample) and B = 2 .74 ± 0.23 (CC sample only) from 
IMorandi etUI (|2007h . 

To compare the normalization in scaling rela- 
tions in the same analytic form and units, we 
convert the SuZIE II normalizations to A = 
A' - log 10 [l g(x,T c )x 3 {e x -l)' 1 ], where 7 = 

2 (^bTcmb) 3 (hc)~ 2 , x is calculated at the SuZIE II ob- 
serving frequency of 145 GHz, and the primes stand for 
the power indices or normalizations from the references 
that w e compare. For normalizations of Morandi et al. 
j2003), A = A' - log 10 (10~ 8 I ) - B' log 10 (7), following 
the same convention as in SuZIE II case. Our values for 
the normalization, A, in both sets are consistent within 
la with the values A = -6.24 ±0.22 from BIMA/OVRO 
(IBonamente et all 120081), -6 .64 < A < -5.82 from 
SuZIE II ([Benson et al.ll200l . and -6.67 < A < -6.14 
for a combined CC+NCC samp le and —6.80 < A < 
-6.35 for a CC-only sample from IMorandi et al.l (l2007V 
Figure [1] shows that the scaling relation based on the 
ASCA/ROSAT parameters has a lower normalization 
than the CTicmdra-based relation due to the systemat- 
ically higher temperatures. The scaling relation is not 
well confined by the ASCA/ROSAT, partly due to larger 
errors and partly due to the fact the scaling is defined 
by only a scatter of the five CC clusters and the single 
NCC cluster A2163. As briefly mentioned in Section[2l if 
A2142 is removed from the Chandra set, since its model 
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TABLE 4 

Scaling relations from X-ray and AMiBA SZE data 



logiofe) 


= A + Blog 10 (x) 




Chandra 






ASCA//J05AT 




X 


y 


A 


B 


xLn(d.O.f-) 


A 


S xL n (d.o.f.) 


T e /keV 
M 25 oo/10 14 M© 
Lx/lO^ergs" 1 


Y 2500 D 2 A E(z)/M P c 2 
Y250oD 2 A E(z)-' 2 / 3 /Mpc 2 
Y 250 oD 2 A E(zf/ 4 /Mpc 2 


°- M4 -0.72 
a o 9 +0.39 

-4 n5+°- 18 


9 90+0.73 

z - z °-0.68 

1 71 +i.oi 

1 -' i -0.64 

0.77±°H 


1.43(4) 
1.82(4) 
6.15(3) 


c q7 +0.67 
°- a ' -0.78 

-5 rn+ - 43 

-4 fiQ+ 28 
4.oy_ 2g 


9 9-, +0.74 
z ' zl -0.64 

1 - au -0.61 
1 11+O.28 
A - L1 -0.29 


3.36(4) 
1.95(4) 
1.69(3) 



Note. — The V2500 — fit uses only five clusters, omi tting A2142 in the Chandra set and omitting A1995 in the ASCA/ROSAT set. 
Xmin gi ves the minimum value of x 2 ■> as defined in Eq. with the corresponding number of degrees of freedom (d.o.f.). Errors are given 

at the 68.3% confidence level. 




k B T e (keV) 

Fig. 1. — The scaling relation between Y2500 and T c . Those 
in solid black show the Chandra-based results. Those in gray 
show the ASCA/ROSAT-based results. Six clusters are la- 
beled as indicated by the legend, with errors represented by 
the boxes. The lines are the best-fit power-law relations. The 
black dashed line is the best fit from Bonamentc ct al. (2008) 
for comparison. 

is somewhat inconsistent with the others, the change on 
the scaling relation is less than 5% because A2142 lies 
close to the best-fit line. 

4.2.2. The Y2500 _ M2500 relation 

The power-law index B based on both sets of X- 
ray parameters are consistent with the self-similar 
model prediction of B = 1.67. Our results also 

agree with th e values of B = 1.66 ± 0.20 from 

BIMA/OVRO (jBonamente et al.ll2008| ). and B = 1.48 ± 
0.39 (CC+NCC sample) an d B = 1.56 ± 0.29 (CC only) 
from yVlorandi et alj (|2007l ) . Our normalization agrees 
with the value A = -5.0 ± 3.0 from BIMA/OVRO 
(jBonamente et al.l2008| ) and is consistent with the ranges 
-5-09 < A < -4.63 (CC+NCC samples) and -5.36 < 
A < -4.91 (CC only) given by IMorandi et all (pOfM . 
Here we convert the normalizations from other studies 
for comparison, such as A = A' + UB' for BIMA/OVRO, 

and A = A' + (B' - 5/3) log 10 E(z ) - log 10 (10- 8 / ) for 
IMorandi et alj (|2007f ). where E(z) is the mean E(z) av- 
eraged over all AMiBA clusters. 

Several analytical and numerical studies demonstrate 
that the integrated SZE signal, I2500 m our case, as a 
measure of the total pressure of inter-clu ster medium is 
an excellent proxy for cluster total mass (Ida Silva et alj 
l200llMotl et alJl2005tlN^^l2TMlHallman et al.ll2007fT 
If this relationship could be measured to high precision 




Fig. 2. — The scaling relation between Y2500 and M25oo- 
Symbols and colors are as define d in Fig. [[] The bl a ck da shed 
line refers to the best-fit from iBon amcntc et alj lj2008f ) for 
comparison. 

at low redshifts, it could then be used to determine the 
masses of high-redshift SZE clusters and test cosmologi- 
cal models. 

4.2.3. The Y2500 ~ Lyi relation 

The power-law index for the Y2500 — Lx relation based 
on the Chandra set (after omitting A2142) is about 1.5cr 
lower than the theoretical value B = 1.25, but is con- 
sistent with the results B = 0.81 ± 0.07 (for CC+NCC 
sample) and B = 0.9 1 ±0.11 (for CC sample) given by 
IMorandi et all (|2007l) . A low power-law index has also 
been observed in numerical si mulations that includ e cool- 
ing or preheating processes (|da Silva et all 12004. Y oc 
Lx). The systematically lower power-law index rela- 
tive to the self-similar model predication seems to im- 
ply that the relation between the SZE signals and X-ray 
luminosities is more sensitive to the radiative content 
outside 100 kpc than other scaling relations. However, 
none of the data points in Fig. [3] lies close to the best- 
fit line and the measure of goodness-of-fit, Xmin/d-O-f-, 
is large (see Tab. 01 . On the other hand, the normal- 
ization agrees with the result —4.35 < A < —4.01 for 
CC+NCC sample within lcr and is bro adly consistent 
with -4.47 < A < -4.36 for CC sample (|Morandi et all 
120071 ) . through the conversion of A = A' - (5/4 - 

B>) log 10 W) ~ logic (I0~ 8 /o) + B>. 

The power-law index and normalization based on the 
ASCA/ROSAT set are consistent with the self-similar 
model. They also a gree with the r e sult b ased on the 
CC sample given by IMorandi et al.1 (|2007l ) within la, 
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Fig. 3. — The scaling relations between Y2500 and Lx- Sym- 
bols and colors are as defined in Fig. [1] The black das hed 
line shows the best-fit relation from M orandi et al.l 1)20071 ) for 
comparison. The luminosities from the ASCA/ROSAT set 
are given without errors. 

but only marginally consistent with those based on the 
CC+NCC sample. We ob served, by co mparing the val- 
ues of different models in lAllenl (|2000l ). that the addi- 
tional component compensating the cooling flow emis- 
sion in Model C will generally reduce the bolometric 
luminosities. If ther e is residual luminous emission, as 
lMorandi~ t aL (200j|) remarked that CC clusters system- 
atically have larger luminosities than NCC ones even if 
the cooling cores have been handled, it would bias high 
the power (slope) and bias low the normalization (inter- 
ception) shown in Fig. [3] in the sense of shifting the CC 
clusters to higher Lx- This is a possible interpretation of 
the discrepancy for the CC+NCC sample, but again the 
scaling relation is not well defined, essentially by a scat- 
ter of CC cluster and a NCC cluster outside the scatter. 

5. DISCUSSIONS 

In the three figures of scaling relations we see that the 
Chandra-based relations are generally better fits than 
the ASCA /ROSAT-based relations, with smaller xLm 
and smaller errors on each data point. Although the 
Y — Lx relation based on the Chandra set has a larger 
scatter, there are no errors available for the luminosi- 
ties of the ASCA/ROSAT set. Among clusters in the 
ASCA/ROSAT set, A2261 and A2390 have X-ray pa- 
rameters of poor quality. A2390 seems to have a biased- 
high gas temperature or a systematically low ^2500- A 
high temperature would lead to a high total mass based 
on the hydrostatic equilibrium equation, and would sim- 
ilarly increase the lumino sity, which is related to the gas 
temperature as Lx oc T e 2 (jMorandi et al.ll2007f ). 

Analytical and numerical studies reveal the fundamen- 
tal incompatibility between /3-model fits to X-ray sur- 
face brig htness profiles and those done with SZE pro - 
files (e.g. IKomatsu fc Seliai3[200lt IHallman et al.ll2007li . 
Both X-ray and SZE fitted model parameters are biased 
due to the isothermal assumption, since the X-ray surface 
brightness and SZE Compton-y parameter have different 
dependences on the cluster temperature profiles. This 
will generate an inconsistency in the model parameters 
based on isothermal /3-model fits. Since observational 
SZE radial profiles are in short supply, X-ray driven pa- 
rameters are often used to constrain the profile shape in 



SZE analysis, consequently leading to a bias in the de- 
rived values of cluster mass or Comptonization parame- 
ter. 

To remedy this problem, wc followed Hallman et al.l 
(|2007t) . Instead of re-fit ting by the un i versal tempera- 
ture profile proposed by Hallm an et al.l (|2007D . we sim- 
ply modify our values of (3, r c and Y2500 by the ratios be- 
tween the values fitted from X-ray data on an isothermal 
/3-model, and the 'true' values obtained from the simula- 
tion. We then re-calculate the scaling relations and these 
corrected results are summarized in Table [SJ It is clear 
that the introduction of correction still keeps the scaling 
relations consistent with the uncorrected results, and the 
previous arguments and discussions are still valid. The 
sc aling relations seem to be insensitive to this correction. 
As IHallman et al.l (|2007D discovered in their study of the 
Y — M gas relation, for example, the correlated changes in 
a /3-model due to the definition of projected radius (^2500 
in our case) tend to cause compensating changes in the 
scaling relation. 

We are aware that the entropy floor present in the cores 
of clusters could give rise to devia t ions fr om self-similar 
scalings (see e.g.. [McCarthy et al.l (|2003allbl )). X-ray ob- 
servations have shown that scaling relations between sev- 
eral cluster observables deviate from the self-similar pre- 
diction, and it has been found that heating and cooling 
act in a similar manner by raising the mean entropy of 
the intracluster gas and, in so me cases, establ i shing a 
core in the entropy profile. In IMcCarthv et al.l ((2003a) 
it was observed that the injection of excess entropy (pre- 
heating) will increase the temperature and reduce the gas 
pressure in the central regions of clusters, especially for 
the low-mass clusters. The scaling relations between the 
central value of the Compton-y parameter, yo, and the 
gas temperature or the total mass are most sensitive to 
the presence of excess entropy, and tend to develop larger 
power indices. Scaling relations involving the integrated 
Compton-y parameter Y show similar behaviors but are 
less sensitive, since the integration tends to reduce the 
effect of the entropy contribution from the cores. 

6. CONCLUSIONS 

In understanding cluster physics and cosmic evolution, 
the study of scaling relations for galaxy clusters is be- 
coming more important since SZE observations such as 
those from the Sunyaev-Zel'dovich Array (SZA) and the 
South Pole Telescope (SPT) will detect clusters which 
are beyond X-ray detection limits, for which SZE scaling 
relations provide the first indications of cluster proper- 
ties. In addition, deviations between the theoretical and 
observational results of the scaling relations can serve to 
examine the non-gravitational processes in the formation 
of clusters, which are not well understood at present. 

As one of the few SZE instruments working at 3-mm, 
the AMiBA experiment observed six Abell clusters dur- 
ing 2007. The derived integrated Compton-y parameters, 
Y250O; are compared to other observations at different 
frequencies, as summarized in Table |31 Our results are 
consistent with those from BIMA/OVRO, but appear to 
show lower Comptonizations than those from SuZIE II. 
We have also investigated the three scaling relations be- 
tween K2500 and the X-ray spectroscopic temperatures, 
total masses within r25oo , and bolometric X-ray luminosi- 
ties. Our results for the scaling relations are summarized 
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TABLE 5 

Scaling relations log 10 (j/) = A + -Blog 10 (x) corrected for isothermal 

/3-MODEL 



X 


y 


A 


£? 


X? nin (d-o.f.) 


T c /keV 
M 2500 /W 14 Mq 
Lx/lO^ergs" 1 


Y 2 moD\E{z) /Mpc 2 
Y 2B00 DlE{z)- 2 / 3 /Mpc 2 
^2500l5iS(^) 9/4 /Mpc 2 


-6 16+ 068 

D - 1D -0.74 

A OO + 0.30 

-4.»i_ .43 

_ 4 24+ - 18 

-0.18 


9 01 +0.75 
z -° -0.70 

u -' a -0.32 


1.74(4) 
1.73(4) 
6.72(3) 



Note. — Errors arc given at the 68.3% confidence level. 



in Table H 

Our power-law indices for the three scaling relations 
are broadly consistent with the self-similar model and 
observational results in the literature, except for that the 
^2500— Lx relation based on C/iandra-derived parameters 
has a slope lower than the expectation of the self-similar 
model, and is sensitive to the chosen set of X-ray parame- 
ters and the treatment of cooling cores. These discrepan- 
cies might indicate either exotic properties for clusters or 
hidden flaws in our SZE data, although the scatter is still 
large, about a factor of two in the integrated Compton- 
y relative to the fit line. The agreement between the 
normalizations found by different workers for our three 
scaling relations seems to support the idea that there is 
no strong scatter in the gas fraction (see Sec. 14. ip . 

In conclusion, the agreement between our results and 
those from the literature provides not only confidence 
for this project but also supports to our understand- 
ing of galaxy clusters. For AMiBA, significant im- 
provements are expected following the expan sion to a 
13-el ement configuration with I.2-m antennas (|Ho et alJ 
2009, AMiBAl3), which will provide better resolution 
and higher sensitivity. The capability of resolving SZE 



clusters will make it possible to mea sure the cluster pro- 
files independent of the X-ray data (jMolnar et al.ll2009T ) 
and to estimate the properties of the clusters which cur- 
rently do not have good X-ray data. 
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